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Pseudopotential study of binding properties of solids within generalized gradient 
approximations: The role of core-valence exchange-correlation. 

M. Fuchs, M. Bockstedte,0 E. Pehlke^l and M. Scheffler 

Fritz- Haber-Institut der Max- Planck- Ges ells chaft, Faradayweg 4-6, D-14195 Berlin- Dahlem, Germany 

(submitted to Phys. Rev. B, July, 11, 1997) 

In ab initio pseudopotential calculations within density-functional theory the nonlinear exchange- 
correlation interaction between valence and core electrons is often treated linearly through the 
pseudopotential. We discuss the accuracy and limitations of this approximation regarding a com- 
parison of the local density approximation (LDA) and generalized gradient approximations (GGA), 
which we find to describe core-valence exchange-correlation markedly different. (1) Evaluating the 
binding properties of a number of typical solids we demonstrate that the pseudopotential approach 
and namely the linearization of core-valence exchange-correlation are both accurate and limited in 
the same way in GGA as in LDA. (2) Examining the practice to carry out GGA calculations us- 
ing pseudopotentials derived within LDA we show that the ensuing results differ significantly from 
those obtained using pseudopotentials derived within GGA. As principal source of these differences 
we identify the distinct behavior of core-valence exchange-correlation in LDA and GGA which, 
accordingly, contributes substantially to the GGA induced changes of calculated binding properties. 
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I. INTRODUCTION 

Generalized gradient approximatioiJS|-| (GGAs) to 
the exchange-coujeJation (XC) energytJ'El in density- 
functional theoryoLl are currently receiving growing in- 
terest as a simple alternative jto^ improve over the local 
density approximation (LDA)l3'tl in ab initio total en- 
ergy calculations. In various respects the GGA proved 
to be more appropriate thaii_the LDA, namefe (1) Bind- 
ing energies of moleculesli30 and solidsESTj turn out 
more accurate, correcting the tendency of the LDA to 
overbinding. (2) Activation energy barriers, e.g. for the 
dissociaiiya^dsorption of H2 on metal and semiconductor 
surfacesLjilZI are in distinctly better accordance with ex- 
periment. Reaction and activation energies for a vaxiely. 
of chemical reactions show a similar improvement .EETlHI 
(3) The relative stability of structural phases seems to be 
predicted mace-cealistically for magnetictJ as well as for 
nonmagneticEj'Ej materials. Bulk structural properties 
are often not improved within GGA. While the lattice 
parameters consistently increase compared to the LDA, 
a closer agreement with experirnenjtaLdata is reported for 
alkali, Sd, and some Ad metalsJj'EirEj An overestimation 
of up to several percent is found however for 5d metals 
and common semiconductors, their bulk modulus^uoiing 
out too small accordingly (typically by < 25%).tZl"E3 

Regarding the understanding of the GGA and further 
advances beyond it, it is important to gain insight into 
the XC related "mechanisms" underlying an eventually 
improved performance, e.g. along the lines pursued in 
Refs. ||,|2|30, and Complementary, careful estimates 
are needed to what degree computational approximations 
in evaluating the total energy interfere with a comparison 
of different XC functionals. With this in mind we exam- 
ine two interrelated issues that have been of persistent 
concern in comparison of the LDA and the GGA based 
on pseudopotential calculations.E^ 

Firstly in how far the nonlinearities associated with 
the XC interaction of core and valence electrons in these 



XC density-functionals can be accounted for by the pseu- 
dopotentials. Computationally it is expedient to treat 
core-valence XC as a part of the pseudopotential and 
thus as if it acted linearly on the (pseudo) valence elec- 
tron density. Within the LDA the transferability of the 
pseudopotentials remains, in most cases, intact under 
this approximation, i.e. good agreement of the results 
of pseudopotential with all-electron calculations may be 
expected without handling these nonlinearities explicitly. 
The extent to which this carries over to GGA's, and thus 
enables a menaningful comparison with the LDA, is un- 
clear at present. Experience with GGA's still needs to be 
built up and previous studies have advanced conflicting 
views on this subject: Examining structural parameters 
of crystalline solids, Juan et aLljo concluded nonlinear 
core corrections for XC were required in pseudopoteru 
tial calculations within the GGA by Perdew and Wangtl 
(PW) even in cases where they are negligible in LDA, like 
bulk Si. On the other hand, Moll et a/.LJ and Dal Corso et 
alkn found the LDA and PW GGA behaved alike in this 
respect. 

Secondly we discuss the role of differences seen in the 
pseudopotentials constructed within the LDA and within 
the GGA. Here we address, on the one hand, in how 
far such differences are small enough to warrant the cir- 
cumventing of a full self-consistent GGA calculation (us- 
ing GGA pseudopotentials) by a computationally simpler 
post-LiDA treatment where the electronic total energy is 
first minimized within LDA (using LDA pseudopoten- 
tials) and then corrected perturbatively for the GGA XC 
energy. On the other hand, the behavior of the pseu- 
dopotentials eventually reflects a different description of 
the core-valence interactions in LDA and GGA. This al- 
lows us to conceive the GGA-related effects separately 
in terms of XC among the valence electrons themselves 
and XC of the valence with the core electrons. Results in 
several works indeed hint that the H)A and GGA might 
differ in this respect: Garcia et alS^ evaluated cohesive 
properties of some metals and semiconductors on the ba- 
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sis of the Becke/Perdew (BP) GGA&Q and a precursor 
to the PW GGA. Dependent on whether the pseudopo- 
tentials screened within the GGA were derived from an 
LDA or GGA calculation of the free atom they obtained, 
in some instances, differing values of the lattice parame- 
ters and cohesive energies. Examining the dissociation of 
silanes using LDA-based pseudopotentials Nachtigall et 
alSB observed noticeably overestimated activation and 
reaction energies compared to the respective all-electron 
approach for various GGAs, but close agreement for the 
LDA. Despite the apparent incongruencies, GGA cak 
tions are still being based on LDA pseudopotentials B3 

To address these issues we investigate the differences 
between the LDA and GGA systematically at each step 
of the pseudopotential approach, the construction of the 
pseudopotentials from atomic calculations and their use 
in polyatomic systems. In turn we evaluate the cohesive 
properties of a set of typical metallic, semiconducting, 
and insulating crystals (Na, Mg, Al, Cu, W, diamond. 
Si, Ge, GaAs, and NaCl), where we apply pseudopoten- 
tials with and without nonlinear core corrections. With 
respect to the role of core- valence XC we establish in how 
far its handling affects the accuracy of pseudopotential 
calculations by comparing our results with available all- 
electron data. We then discuss the related need for the 
consistent use of the same XC scheme at all points of a 
pseudopotential calculation and comment on the contri- 
bution to the GGA induced changes of LDA results for 
cohesive properties driven by differences between LDA 
and GGA core- valence XC. 

Concerning proposals for GGAs we present results for 
the PW and the earlier BP schemes. Both are variants 
of the generic type 
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depending locally on the electronic density n{r) and its 
gradient and yielding a local XC potential Vxc{^) — 
SExc[n]/dn{r) like in case of the LDA. These schemes 
are widely used in present day applications, and remain 
of interest as a starting point in recent nonlocal hybrid 
XC schemes-expected to further improve over GGA type 
functionals.Ea Accurate all-electron results are available 
for those GGA schemes and serve as a rigorous reference 
for the pseudopotential calculations in this study. The 
PW GGA is derived basically from first principles, com- 
bining the gradient expansions of the exchange and corre- 
lation holes of a non-uniform electron gas with real-space 
truncations to enforce constraints imposed by properties 
of the physical XC hole. While the BP GGA may be 
deemed to be somewhat more heuristic as it relies also 
on fitted parameters, it has been yielding results close 
to those of the PW GGA, at least in all-electron cal- 
culations. In addition we have considered the recently 
proposed GGA by Perdew, Burke and Ernzerhof (PBE) 
which is regarded to be conceptually more concise than 
the PW GGA but is expected to perform essentially sim- 
ilarlyEj In the pseudopotential calculations for the prop- 
erties addressed here we have found the PiB|E and the 
PW GGA to yield nearly equivalent results Jlj hence our 
conclusions for the PW GGA hold for the PBE GGA as 
well. 



The remainder of this paper is organized as follows: In 
Sec. H we briefly review and discuss the relevant formal 
aspects of pseudopotential calculations. Technical char- 
acteristics of our calculations are outlined in Sec. |lll|. In 
Sec. 
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' we report our results, and put them in perspective 
with our considerations from Sec. ||. Section ^ summa- 
rizes our conlcusions. Atomic units are used throughout 
unless indicated otherwise. 



II. GENERAL CONSIDERATIONS 

In pseudopotential calculations the total energy is for- 
mally treated as a functional of the valence charge den- 
sity alone, with the pseudopotentials accounting for the 
interaction of the valence electrons with the nuclei and 
with the core electrons - namely for Pauli repulsion, elec- 
trostatic and XjQiateractions - to within the frozen core 
approximation.[i2lc3 Substituting the GGA for the LDA 
modifies the treatment not only of the XC interactions 
of the valence electrons among themselves but also that 
of the core-valence (CV) interactions. In order to treat 
all interactions within one and the same XC scheme the 
pseudopotentials to be employed in a GGA calculation in 
principle ought to be generated consistently within the 
same GGA as well, rather than within, say, the LDA. 
In the following we discuss the relevance of this "pseu- 
dopotential consistency" to total energy calculations in 
LDA and GGA. Within the pseudopotential framework 
the GGA total energy functional reads 

occ 
i 

(2) 

where the various terms denote the noninteracting ki- 
netic energy, the Hartree energy, the XC energy and the 
potential energy of the valence electrons, represented by 
the pseudo wavefunctions '0i(r) and the corresponding 
charge density n(r) — ^2°'^'^ |'0i(r)p in the presence of the 
ion cores, represented by GGA pseudopotentials, V^ciga^ 
The LDA counterpart to (||) is obtained by substitut- 
ing the XC energy -E^xc^ ^^'^ LDA pseudopotentials 
prLDA j-Qj. ^]-^g respective GGA entities. 

Now the groundstate energies in GGA and LDA can be 
readily compared with the help of a perturbative analysis 
of the total energy functionals E^f^ [n] and E^^^ [n] , at 
any given set of ionic positions. Around the respective 
stationary groundstates, characterized by the densities 
^GGA n^^^^ the variational principle implies that 
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and likewise in LDA. Supposed GGA and LDA yield sim- 
ilar densities the difference of their groundstate energies, 

El'^^[n^^^], can be expressed 



SEtot — j^tot 
by virtue of (t 
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i.e. simply in terms of the density and wavefunctions ob- 
tained within the LDA, using the LDA pseudopotential, 
yLDA Accordingly the GGA modifies the (pseudo) to- 
tal energy in two ways: (i) By the direct difference of 
the XC energies SExc = E^§^[n^^^] - E^^^[n^^^]. 
This term corresponds to the often applied a posteriori 
gradient-correction scheme where, at given ionic coordi- 
nates, the density is evaluated self-consistently from the 
LDA XC potential and then used to compute the total en- 
ergy with the GGA XC energy functional. For the GGAs 
considered the SExc is negative, and vanishes in the lim- 
iting case of the homogeneous electron gas. Typically 
its magnitude increases with the degree of inhomogene- 
ity of the system at hand,c3 the GGA correction to the 
LDA XC energy being larger for a free atom or molecule 
than in a solid, (ii) By the potential energy correction, 
SEv = - V^^^l^jjf^^), which arises 

as a consequence of the pseudopotential approximation 
and eventually reflects the differences in the behavior of 
the CV interactions in LDA and GGA. Note that in an 
all-electron formulation this term would be absent alto- 
gether. Clearly the potential energy correction is missed 
when LDA pseudopotentials are carried over to GGA cal- 
culations, giving rise to a "portability" error compared to 
the consistent GGA calculation using GGA pseudopoten- 
tials. We shall demonstrate in Section IV that SEy does 



not cancel out when total energy differences are consid- 
ered but is in general of similar importance as SExc for 
quantitative tests of the GGA within the pseudopotential 
framework. 

For an understanding of the differences between the 
GGA and LDA it is worthwile to examine more closely 
the various contributions to the CV interactions the pseu- 
dopotentials mediate. In the following we identify aBd, 
discuss these for norm-conserving, pseudopotentials ,e2I 
constructed by standard schemesEJ'El from atomic all- 
electron calculations. As a canonical first step these algo- 
rithms generate angular momentum dependent screened 
pseudopotentials Vj'^'^fno] from a particular reference con- 
figuration, e.g. the groundstate of the neutral atom, as- 
suming spherical screening. These act as effective poten- 
tials on the atomic pseudo valence states via the radial 
Schrodinger equations, 
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+ Vf''[no;r]-ei rRi{r)^Q. (5) 



The Vf [hq] contain a common spherical screening po- 
tential which is self-consistent with the total atomic 
charge density no(r), comprised of the (pseudo) valence 
density Uq and the core charge density Uq obtained from 
the all-electron core states. The effective potentials can 
be decomposed rigorously into the Hartree potentials Vh 
and the XC potential due to the valence and core elec- 
trons and an angular momentum dependent bare poten- 
tial Vi^'^'^^ which conveys the nuclear attraction and the 
Pauli repulsion due to the core states; for an arbitrary 
valence configuration one has 

y,^*^[n;r] = y^"(r) + yHK;r] 

+14jK;r] + l/xcK+nS;r], (6) 



which, in the reference configuration [n = tIq + Ug), 
reduces of course to the screened pseudopotentials. 
Through the nonlinearity of the XC potential in the den- 
sity the effective potential retains a dependence on the 
total rather than on the valence density alone as it is 
ultimately a prerequisite for an efficient plane wave rep- 
resentation. Customarily a further separation in terms of 
frozen core and variable valence contributions is accom- 
plished by "linearizing" the XC interaction taking 

+ {Vxc <,r]~ Vxc [n^o+ «S ; ^] ) , (7) 
where the partial core density nQ{r) serves as a control 



parameter. Choosing n"" 



the screened pseudopo- 



tentials, and thus the atomic properties in the reference 
configuration, are correctly recovered. Now those terms 
on the RHS of ^ which are independent of the valence 
density define the usual pseudopotentials that are to be 
transferred to and screened according to the environment 
of one's target system. Applying they read 



(^xc K+ no;r]- Vxc "S; ^]) 



(8) 



The last term here comprises the core- valence XC inter- 
action 

AVxcir) - ("^5ccK+ "S;^] - VxcK+ (9) 

as represented by the ionic pseudopotential. Below we 
demonstrate that it is the key quantity to understand the 
differences between the ionic pseudopotentials in LDA 
and GGA which, in practice, are defined simply by "un- 
screening" the screened potentials according to 

V,{r) = V'''[no;r] - Vii[n-„;r] ~ FxcK+^C;^]- (10) 

with all quantities evaluated within the respective XC 
scheme. Note that the transformation Eq. (0) turns the 
core-valence XC energy into a linear functional of the 
valence density that is absorbed in the pseudopotential 
contribution to the total energy instead of being treated 
as a part of the XC energy itself. By experience the com- 
plete core-valence linearization, ng(r) = 0, has proven 
to be accurate for the majority of applications within 
the LDA. It is expected to be justified for local func- 
tionals like the LDA and also the GGA if the overlap of 
the core and valence charge densities does not substan- 
tially change whenever chemical bonds are formed or al- 
tered. Formally the nonlinear core-valence XC could be 
regarded exactly, taking ng(r) = ng(r), and, correspond- 
ingly, Exc =■ Exc[n + TIq] in the total energy functional 
(Q), where rig denotes the core charge density as com- 
pounded from the frozen atomic core charge densities. 
Since the core states are strongly localized and sharply 
peaked such a choice is beyond the realm of a plane wave 
representation however. If a complete linearization of CV 
XCl4)roves insufficient, e.g. in calculations of alkali met- 
alsLJ or involving spin-polarization,Lj the nonlinearities 
can still be captured adequately in the chemically most 
important interatomic regions with the help of af^artial 
core density, as was first realized by Louie et a^J£3 It is 
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tailored to coincide with the full core charge density be- 
yond a suitable cutoff radius Tj. but avoids the sharply 
peaked structure close to the nucleus by a smooth cutoff 
function b{r) chosen largely at computational expediency, 
cf. Rcf. 



and Sec. 

ng(r) 



IV, 



for r > Tr 



6(r)n3(r) for r < Vc, with b{r) < I. 



(11) 



Note that the CV XC component of the ionic pseudopo- 
tentials vanishes beyond rc- 

It is well understood that the GGA does not substan- 
tially alter the wavefunctions and the spectrum of atomic 
valence states compared_to the LDA, these are bound too 
weakly in both schemes,t2l mainly because the XC poten- 
tials of these schemes insufficiently cancel the repulsive 
contribution from the electrons' self- interaction to the 
Hartree potential. Consequently the effective potentials 
for the pseudo valence states ought to be close for 
both XC schemes. Indeed the screened LDA and GGA 
pseudopotentials are barely distinguishable by a simple 
visual inspection as can be seen, e.g., for germanium and 
copper in Fig. |^. 
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FIG. 1. Screened pseudopotentials within BP GGA and 
LDA for germanium and copper. On the scale of these plots 
the GGA and LDA pseudopotentials lie one on top of each 
other. Shown are Troullier-Martins pseudopotentials with 
cutoff radii rs,p — 1.9 bohr and Vd ~ 2.3 bohr (Ge), and 
rs,d = 2.0 bohr and = 2.3 bohr (Cu). 

Turning to the unscreened, ionic pseudopotentials ac- 
tually used for calculations more pronounced deviations 
between LDA and GGA pseudopotentials emerge. These 
may by easily analyzed in terms of the various pseudopo- 
tential components given by Eqs. §) and @) which con- 
tribute to the difference 5Viir) = V;^G^(r) - V^^^^ir), 
where the superscripts indicate the type of XC employed 
in constructing the pseudopotentials. This decomposi- 
tion is illustrated in Fig. |[ weighting all differences with 
the r-dependent volume element. To highlight the role of 
the individual contributions we distinguish three cases: 

(i) completely linearized XC, n.o(^) = 0- 

(ii) approximate account of nonlinear CV XC, employ- 
ing a partial core charge density identical with the 



full one outside Tc = 1.3 bohr. 

(iii) full account of nonlinear CV XC, taking 7T.g(r) = 
n[j(r). 

Case (iii) serves to identify the genuine difference poten- 
tial, due to the unlike bare and core-valence Hartree po- 
tentials. These originate just from the small differences 
of the self-consistent atomic orbitals in LDA and GGA 
and thus remain the same in (i) and (ii) where CV XC is 
approximated. The core-valence Hartree potentials are 
seen to make a small positive contribution, 

SVn{r) = l%K(GGA);r] - HjK(LDA);r]. 
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FIG. 2. Difference of the ionic pseudopotentials in BP 
GGA and LDA for germanium and copper, corresponding 
to the screened potentials given in Fig. ||. Shown are the 
XC, Hartree, and bare potential contributions, cf. Eq. (^, as 
discussed in the text. Both panels use the same scale for the 
ordinate. The cutoff radius of the partial core charge densities 
was chosen as 1.3 bohr. 

This reflects the fact that the core states are more 
tightly bound in GGA than in LDA, some.\V|lji at enhanc- 
ing the electrostatic screening of the nucleLcacj Still SVu 
is found to be weaker and faster decaying than the (os- 
cillatory) differences in the bare potentials. In case (i) 
the genuine difference potential is superimposed with a 
long ranged, repulsive hump. This feature clearly signi- 
fies the distinct analytical behavior of the CV XC poten- 
tials, Eq. (H), in LDA and GGA, rather than differences 
in the self-consistent charge densities. Notably for germa- 
nium it attains its maximum around and stretches well 
beyond the maximum of the valence charge density up 
to radii that correspond to mid-bond positions. Includ- 
ing a partial core charge density, case (ii), eliminates by 
construction the core-valence XC potentials outside the 
respective cutoff radius in both LDA and GGA. In this 
region one therefore recovers the more short ranged, and 
in case of germanium weaker, genuine difference poten- 
tial. Inside there remains some interference of genuine 
and core-valence XC related differences as partial and 
full core charge density deviate from each other. 

The above discussion suggests that the distinct core- 
valence XC interaction in LDA and GGA is a prime 
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source of the differences between the pseudopotentials. 
As the LDA pseudopotentials and their GGA counter- 
parts differ even in interatomic regions they are to be ex- 
pected to perform unhke in a pseudopotential calculation 
which employs the GGA for the XC energy, germanium 
with linearized CV XC being a generic example. By ex- 
plicitly considering nonlinear CV XC, or, as exemplified 
by copper, by including more semicore states as valence 
states the difference in CV XC is removed from the pseu- 
dopotentials and instead taken into account through the 
XC energy functionals themselves. In this case LDA and 
GGA pseudopotentials should thus behave more alike, 
provided of course the cooperative genuine differences are 
negligibly small themselves. In Sec. IV we substantiate 
these aspects quantitatively. 



III. COMPUTATIONAL METHOD 

From a practitioner's point of view, GGAs are readily 
incorporated in plane wave based schemes: the deriva- 
tives of the density needed to compute XC energy and 
potential in position space are evaluated from the recip- 
rocal space representation of the density and transformed 
to position space using Fourier transformations, the con- 
vergence of all relevant quantities being controlled - like 
in case of LDA - through the plane wave basis size. The 
construction of norm-conserving pseudopotentials within 
GGA's proceeds entirely parallel to the one in LDA. The 
necessary radial density gradients may be inferred, e.g., 
directly from the derivatives of the radial wapifunctions. 

We have constructed the pseudopotentialsc£l based on 
a scalar-relativistic atoirvie, calculation using the scheme 
of TrouUier and Martins .eJ Core and valence states were 
partitioned as usual, i.e. retaining only the uppermost 
occupied s and p states as valences, except for Cu and 
W where the 3d(5d) states need to be included in the 
valence space. The resulting semilocal potentials were 
further transformed into fully separable representations 
of the Kleinman-Bylander kindE3 In case of nonlinear CV 
XC we used a cuspless polynomial to represent the partial 
core charge density inside the cutoff radius. Continuity 
of the density up to its third derivative is enforced to en- 
sure that the GGA XC potential joins smoothly. Various 
tests, carried out for the free pseudo atoms and described 
in the appendix, indicate comparable transferability for 
the GGA and LDA pseudopotentials. The ionic pseu- 
dopotentials are tabulated and transferred without any 
intermediate fitting to the plane \w/e calculation. We 
have refrained from any smoothingp3 of the ionic GGA 
pseudopotentials which on occasion display short-ranged 
oscillations, mostly for the PW GGA.EjS These corre- 
spond to a plane wave energy regime where the kinetic 
energy is dominating all other energy contributions and 
may, therefore, be conceived to be physically negligible. 
Care is required though to leave the relevant low Fourier 
components intact if smoothing is performed. With a 
numerical tabulation this is attained in an unbiased, sys- 
tematic manner through the basis size cutoff. 

We have computedE^ the total energy per atom in the 
bulk systems varying the lattice constant within about 



±5% of the respective equilibrium value. |-Jitting these 
energies to Murnaghan's equation of stateES we obtained 
the equilibrium values of the lattice parameters and the 
total energy per atom. The cohesive energy was deter- 
mined by subtracting the latter from the total energy 
of the spin-saturated spherical (pseudo) atom. To cor- 
rect this value for neglected contributions due to the 
spin-polarization of the atomic groundstate, we added 
the difference of the total energies of the spin-polarized 
and spin-saturated all-electron atom within the respec- 
tive XC scheme. Corrections of the theoretical values of 
the cohesive energy for the phonon zero-point energies 
are disregarded, they amount to «180 meV for diamond 
and are, expected to stay below ss60 meV for the other 
solids.El The Brillouin zone sampling for the bulk sys- 
tems was carried out using 6x6x6 (diamond, NaCl), 
8x8x8 (Al, Si, Ge, GaAs)pand 10 x 10 x 10 (Na, Cu, W) 
meshes of special k-points.o For evaluating the cohesive 
energies we chose a plane wave cutoff energy of 50 Ry for 
all crystals other than diamond, Cu, and W for which 
we used 100 Ry. The respective structural parameters 
were determined with roughly two thirds of these values. 
These computational parameters allow for a numerical 
precision of better than 0.5% for the lattice coBstants 
and better than 50 meV for the binding energies 



IV. RESULTS AND DISCUSSION 

Two sets of GGA calculations were performed where 
we adopted either the consistent GGA approach, employ- 
ing the GGA for both plane wave calculation and con- 
struction of the the pseudopotentials, or, by contrast, the 
inconsistent GGA approach, employing the GGA for the 
plane wave calculation but LDA for the construction of 
the the pseudopotentials. The results of our calculations 



are compiled along with reference data in the Tables IV 



VIII . By a comparison with all-electron data and calcula- 



tions including nonlinear CV XC we first demonstrate the 
pseudopotentials to be as transferable within the GGA 
as in the LDA. In particular transferability is not more 
stringently limited by nonlinear CV XC in GGA than in 
LDA. This provides the frame of reference for our sub- 
sequent examination of the consistent and inconsistent 
GGA approaches which are found to be inequivalent in- 
deed. Following our discussion of Sec. g we then show 
that the perturbative potential energy correction yields 
a good quantitative account of the differences between 
these approaches. Together with the fact that the dis- 
crepancies are essentially eliminated once nonlinear CV 
XC is included, this allows us to identify these discrepan- 
cies as a manifestation of the distinct behavior of CV XC 
in LDA and GGA. Namely we find CV XC to contribute 
significantly to the correction of the binding energies and 
lattice parameters induced by the GGA. 



A. Pseudopotential transferability within GGA 

Possible uncertainties rooted in the pseudopotential 
approximation itself are properly distinguished from ef- 
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fects due to the use of different XC functionals by a com- 
parison with aU-electron data. To this extent we show 
in Fig. ^ the relative error of the lattice constant with 
respect to its experimental value based on a co mpil ation 
of results from recent all-electron calculationsEaES and 
as obtained in the present pseudopotential framework, 
accounting for nonlinear CV XC and working with the 
consistent approach to PW GGA. It can be seen that in 
LDA as well as in PW GGA the results of the pseudopo- 
tential and the all-electron method agree on the order of 
or better than 1%. In particular both methods yield vir- 
tually the same lattice expansion due to the GGA com- 
pared to LDA. Merely in case of Al the agreement is not 
fully quantitative. Concerning the bulk moduli we note 
that the pseudopotential calculations reproduce the re- 
duction of the LDA values due to the GGA as obtained 
from the all-electron calculations. The residual discrep- 
ancies between pseudopotential and all-electron results 
for the bulk moduli of Ge and GaAs in LDA as well as 
in GGA may be seen as a reminder that a more accurate 
treatment requires the extended 3|£Lfitates of Ga and Ge 
to be considered as valence states.Ea 
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FIG. 3. Top panel: Relative error of the bulk lattice con- 
stant evaluated within LDA and GGA with respect to the 
experimental value. Filled squares refer to the present work, 
open circles to all-electron results of Refs. and ^ Bottom 
panel: Increase of the bulk lattice constant (Sao) from the 
LDA to the GGA value relative to the experimental value for 
either method. 

The routine disregard of nonlinear CV XC entails no 
alterations of the calculated lattice properties compared 
to the results obtained with explicit account of nonlin- 
ear CV XC for diamond, Al, and Si. For Ge and GaAs 
the lattice constants are reduced by wl% in LDA and 
changed somewhat less in GGA. The values of the lattice 
constants of sodium metal turn out smaller by ~2% in 
LDA and change again less for the GGAs. Those of NaCl 
agree with the experimental values at the usual level, but 
turn out unrealistic with linear CV XC. It is well estab- 
lished that an explicit account of nonlinear CV XC is es- 
sential in order to predict dependable lattice properties 
for compounds of alkali metals within the LDA.L3 Our 
results for NaCl suggest that this conclusion applies to 
the GGAs as well. Such a behavior seems reasonable, as 



in view of the relatively easily polarizable valence shell of 
Na its core- valence overlap in metallic sodium or NaCl is 
likely to depart considerably from the one in the isolated 
Na atom so that an explicit account of the ensuing nonlin- 
ear changes of core- valence XC becomes indispensable at 
large. The presence of slight differences in the calculated 
lattice properties with and without explicit nonlinear CV 
XC for GaAs, Ge, and W can be similarly conceived as a 
signature of the extended core charge densities in these 
atoms compared to C, Al and Si, for which such dif- 
ferences are not observed. They are likewise absent in 
copper, where the semicore 3 rf-electrons are considered 
as valence states so that their XC interactions with the 
4s-electrons are incorporated exactly. Thus we altogether 
find the neglect of nonlinear CV XC within the GGAs of 
similar importance to the transferability of the pseudopo- 
tentials as within LDA. In all cases where nonlinear CV 
XC is dispensable in LDA it proved to be negligible in 
the GGAs as well. Judged by the systems and properties 
considered we clearly find the pseudopotential approach 
itself to be equally applicable and accurate in GGA as in 
LDA. 

Our results compare well with those of previous pseu- 
dopotential calculations reported in Refs. ^ (BP GGA) 
and |7| (BP and PW GGA), and Ref. g (PW GGA). The 
severe overcorrection of the lattice constants of Al, Si, Ge, 
and GaAs in case of PW GGA without nonlinear CV XC 
found by Juan and KaxirasH is not confirmed here. Sim- 
ilar to Ref. ^ our findings do namely not support the 
conjecture put forward by these authors in Ref. |3j that 
an explicit treatment of nonlinear CV XC were necessary 
in order to arrive at transferable pseudopotentials within 
the PW GGA. 

For the present sytems we find both GGAs to predict 
closely agreeing structural properties with only immate- 
rial differences. Compared to the LDA we note enhanced 
agreement with experimental data regarding the lattice 
constants of Na, Mg, Al, and Si. For the other materials 
the GGA functionals overestimate the lattice constants to 
a similar degree as LDA underestimates them. The bulk 
moduli within the GGAs are predicted in good accor- 
dance with experiment only for Na and Cu, for the other 
materials they are clearly underestimated, in particular 
(by up to w 25%) for the semiconductors. In all cases 
the PW GGA yields cohesive energies in close agreement 
with experimental figures and corrects the overbinding 
of the LDA. Observe however that in GGA the values 
of the atomic energy and hence the cohesive energy were 
still lowered by up to several tenth of an jSiZ—by allow- 
ing for nonspherical groundstate densities J1j'E3 possibly 
unveilling a slight "underbinding" indicated already by 
the overestimate of the lattice parameters. We note the 
BP GGA to yield cohesive energies systematically lower 
than in the PW GGA suggesting a slightly weaker bind- 
ing than the PW GGA. 



B. LDA vs. GGA pseudopotentials in GGA 
calculations 

Having reassured ourselves of the validity of the pseu- 
dopotential ansatz itself we now turn to an account of 
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the inconsistent GGA approach where LDA rather than 
GGA pseudopotentials are employed. As argued in Sec- 
tion H, LDA and GGA pseudopotentials exhibit substan- 
tial differences even in the interatomic regions of molec- 
ular or crystalline compounds. This eventually implies 
that LDA and GGA pseudopotentials perform differently 
when they are combined with the GGA XC energy func- 
tional. 

It is evident from the results listed in Tables | - VIII 
that the outcome of the inconsistent approach depends 
sensitively on the handling of core- valence XC. We shall 
address linear CV XC first. Within it, the inconsistent 
GGA approach yields a description of the cohesive prop- 
erties clearly disparate to that obtained from the con- 
sistent one: the characteristic lattice expansion in either 
BP or PW GGA does not occur, and instead the lat- 
tice parameters closely resemble their LDA values. Like- 
wise the cohesive energies turn out larger by 0.3 eV to 
0.7 eV, amounting to «10% (diamond) up to «50% (Ge) 
of the GGA correction to the cohesive energy found with 
GGA PPs. A somewhat modified behavior is observed 
for copper where the lattice parameters are close to those 
from the consistent GGA approach and the cohesive en- 
ergy correction turns out larger by about 0.2 eV than 
with GGA pseudopotentials. We find the consistent and 
the inconsistent GGA approach to both yield congruous 
descriptions of the binding properties once the calcula- 
tions include nonlinear CV XC: the inconsistent GGA 
approach uniformly recovers the typical lattice expansion 
as well as the decrease of the cohesive energy. The only 
incongruencies found between the two approaches con- 
cern residual deviations of about 0.2 eV for the cohesive 
energies of diamond and copper. 

To further discuss the GGA correction of the cohe- 
sive energy, SEi,, dependent on the choices for pseudopo- 
tential and the treatment of CV XC, we examine the 
constituent corrections of the total energies for (pseudo) 
atom and solid separately by the decomposition 



E; 



LDA 



(12) 



SEl 



alid 



making use of the the perturbative analysis of Section ||. 
Following Eq. (^ the GGA entails a twofold change of the 
total energy compared to the LDA: the direct correction 
of the XC energy, SExc , and the potential energy correc- 
tion, SE^. Now the consistent approach comprises both 
corrections so that the change in total energy is given 
by SEtot — SExc + SEy. By contrast the inconsistent 
approach neglects the difference of the LDA and GGA 
pseudopotentials so that SEy vanishes and the change in 



total energy is limited to SEtot — SExc ■ In Table Xl 



detail the various terms for some exemplary cases. 

We first of all see the perturbative treatment to be well 
justified as it closely reproduces the total energy correc- 
tions extracted from the respective self-consistent LDA 
and GGA calculations in both the atoms and the solids to 
within 0.02 eV. Adopting linear CV XC, the XC energy 
correction in each case is found to account only partly for 
the GGA induced change of the cohesive energy. Instead 
a substantial fraction must be attributed to the potential 
energy correction and is thus indeed effected by the dif- 



ference of the LDA and GGA pseudopotentials. Acting 
in a like manner as SExc^ SEy results in a decrease of 
the cohesive energy, except for copper where we observe 
the above mentioned slight enlargement. Taking non- 
linear CV XC into consideration, the magnitude of the 
potential energy correction is reduced compared to linear 
CV XC. Importantly, 6Ey no longer contributes to the 
change of the binding energy which is instead captured 
completely by the XC energy correction for Na, Si and 
Ge, shown in detail in Table XL It does retain signifi- 



cance in case of diamond and copper however. Save for 
these, the inconsistent and consistent approaches with 
explicit nonlinear CV XC thus produce alike values for 
the binding energies. 

We note the lattice expansion observed upon switch- 
ing over from LDA to GGA to be consistent with the re- 
pulsive character of the difference potential between the 
LDA and GGA pseudopotentials seen in the real space 
inspection of the respective pseudopotentials, cf. Fig. ||. 
Similarly the potential energy correction turns out posi- 
tive, and more so in the solid where valence charge accu- 
mulates in the bonding region. Once nonlinear CV XC 
is taken into account, the XC related differences outside 
the core region are eliminated. What remains are es- 
sentially the genuine differences of the LDA and GGA 
pseudopotentials, reflecting the different description of 
the core states in the respective XC schemes. In prin- 
ciple the more tightly bound core states in the GGA 
should make the Pauli and the Coulomb repulsion more 
short-ranged compared to the LDA, adding some repul- 
sion about the ion sites and some attraction at interme- 
diate distances. The details are certainly quite species 
dependent and moreover one cannot rule out some inter- 
ference from residual XC related differences, e.g. due to 
the unlike partial and full core densities, frustrating any 
a priori estimate of their contribution to the potential 
energy correction SEy. Nevertheless as the genuine dif- 
ferences are small and confined to the immediate vicinity 
of ion sites they should be rather inconsequential to the 
calculation of total energy differences. Such a scenario 
is conceivable and consistent with our results for Na, Al, 
Si, etc., but has its limitations as major portions of the 
valence charge density reside and adjust to charge trans- 
fer close to the ion sites. This clearly applies for the 2p 
states of the first row elements like C and the 3d transi- 
tion metals like Cu, which take on their maxima in the 
domain of the genuine differences. Hence we find the 
potential energy correction for these elements to be only 
partly conditioned by the treatment of nonlinear CV XC. 
In case of Cu the 3d electrons are considered as valence 
rather than as core electrons so that their XC interac- 
tions with the 4s electrons are accounted for exactly and 
not linearly through pseudopotentials. For such a core- 
valence partitioning it is readily verified by inspection, cf. 
Fig. 1^, that the differences of LDA and GGA pseudopo- 
tentials are primarily of the genuine kind and thus quite 
independent of a further account of nonlinear CV XC 
with still deeper core states. Accordingly we obtain an 
actually equivalent description of the bulk lattice param- 
eters in GGA with either LDA or GGA pseudopotentials, 
with the potential energy correction to the cohesive en- 
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ergy being of the order of 0.2 eV within both Hnear and 
nonhnear CV XC. We would Hke to point to an analogous 
observation made in an all-electron method using pseu- 
dopotentials, where the full core density is retained: in- 
vestigating Fe within the PW GGA, Cho et alE3 reported 
nearly identical structural and magnetic parameters with 
either LDA or GGA based pseudopotentials. 

V. SUMMARY AND CONCLUSIONS 

In conclusion we have shown that GGA pseudopoten- 
tials indeed convey a substantial share of the GGA's cor- 
rections over LDA. Accordingly we deem the consistent 
use of the GGA in the application of the pseudopotential 
and their construction to be generally a key requirement 
to attain GGA quantities equivalent to those obtained 
within GGA all-electron methods. By contrast, inequiv- 
alent results arise when the XC energy is treated in GGA 
but the pseudopotentials are taken, inconsistently, to be 
the same ones as in LDA. We have found such "portabil- 
ity errors" to be most significant when the XC interaction 
of core and electrons is treated linearly as a component of 
the pseudopotential, but less important when the nonlin- 
ear core-valence XC interaction is incorporated properly 
into the XC energy functional of the valence electrons it- 
self employing a partial core density. The precise agree- 
ment of the results for the cohesive energies from our per- 
turbative and self-co]aS|istent calculations conforms with 
the the common loretj'cEl that self-consistency has only 
a small effect on the value of GGA total energies and 
differences thereof. Instead these can be accurately eval- 
uated with the LDA wavefunctions and charge density. 
Our analysis shows however that such an a posteriori 
GGA scheme within the pseudopotential framework must 
treat the difference in the valence-valence XC energies 
and the ionic pseudopotentials on an equal footing, as 
given in Eq. (Q): either carrying out the initial LDA cal- 
culation with LDA pseudopotentials and adding to the 
total energy the correction SExc + SEy or, equivalently, 
doing the LDA calculation with GGA pseudopotentials 
{6Ey = 0) and adding just 6Exc to the total energy. 

As the differences of the pseudopotentials originate 
from the distinct core- valence XC potentials in LDA and 
GGA we moreover understand our findings as evidence 
that the bond softening in GGA is directly related to 
a stronger XC repulsion between the valence and up- 
per core states than in the LDA. The reduction of the 
binding energy by the GGA on the other hand appears 
to a larger extent due to describing the XC of the va- 
lence electrons among themselves within GGA instead 
of LDA. The notion of a more repulsive nature of GGA 
core-valence XC agrees with and qualifies earlier obser- 
vations that the GGA corrections to bonding properties 
in solids arise mainly from the immediate vicpity of the 
ions rather than from the interstitial regions.EZl Likewise 
it is supported by the fact that the GGA XC potential of 
atoms like the exact XC potential is superimposed with 
a peaked structure which acts repulsively rat shell bound- 
aries compared to the LDA XC potentiaLE^I Interestingly 
our findings suggest that inadequacies in the description 



of core-valence XC are an important aspect of deficien- 
cies in the description of chemical bonds within either 
LDA or GGA. 

In concluding we note that a conceptual parallel of 
the ("inconsistent") combination of LDA pseudopoten- 
tials with GGA XC is encountered in wavefunction- 
based many-body methods^ such as quantum Monte- 
Carlo (QMC) simulationsES In applications QMC has 
been mostly combined with pseuikipotentials derived 
from effective one-particle schemes.L^I Thereby the inter- 
actions among the valence electrons are described exactly 
whereas the effects of the core electrons are dealt with in 
an approximate manner, say, on the level of LDA. In prin- 
ciple such QMC calculations provide an exact reference 
against which we could check the performance of approx- 
imate XC schemes like the GGA for the valence electrons 
alone. Indeed a survey of the lite£a|t«re indicates that co- 
hesive energies for diamond, Si,E30 and GeE3 obtained 
with QMC and LDA pseudopotentials are in significantly 
closer agreement with experimental figures than are our 
GGA results using LDA pseudopotentials. At least for 
these cases this raises the question whether the GGA 
affords a better description of XC among all electrons, 
then yielding highly accurate binding energies, than of 
XC among the valence electrons alone. On the other 
hand, one is well aware that the use of LDA pseudopoten- 
tials in QMC simulations introduces some uncertainty io. 
the QMC values for the cohesive energy themselves J£3'E3 
quite analogous to the "portability errors" we have found 
in GGA calculations using LDA pseudopotentials to be 
of the order of up to some tenths of an eV. On this level 
of accuracy it is then clearly desirable to obtain accu- 
rate estimates of the systematic uncertainties related to 
the use of pseudopotentials derived from DFT in exact 
methods like e.g. QMC simulations as well in order to 
facilitate a quantitative assessment of approximate XC 
schemes like the GGA. 



VI. APPENDIX 

In this appendix we present some tests on the trans- 
ferability of our GGA pseudopotentials compared to the 
LDA ones. These serve to further corroborate that GGA 
and LDA pseudopotentials show a similar inherent trans- 
ferability but exhibit significant differences due to core- 
valence (CV) XC, as discussed in Sec. ||. As is rather well 
established transferable pseudopotentials should closely 
preserve: (1) The all-electron atomic scattering proper- 
ties as given by the logarithmic derivatives at some radius 
outside the core region over the range of valence energies 
relevant to chemical bonding, say up to ±1 hartree about 
the rpfprpnrp energies. (2) The all-electron atomic hard- 
nesSjHEj i.e. reproduce total energy and eigenvalues for 
excited atomic configurations, to within the accuracy of 
the the underlying frozen-core approximation. 

In Fig. ^ we show the logarithmic derivatives evalu- 
ated with screened pseudopotentials, cf. Eq. (|5|), taking 
germanium as an example. Good agreement with the 
respective all-electron logarithmic derivatives, to be ex- 
pected from the norm-conservation constraints, is con- 
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firmed for both LDA as well as GGA pseudopotentials 
in the serailocal and also in the Kleinman-Bylander rep- 
resentation. For the latter we have additionally }ie|rified 
the absence of ghost states following Gonze et alSB 
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FIG. 4. Logarithmic derivatives R[{e)/Ri{e) vs. energy 
e for the germanium atom (at r — 2.4 bohr, and with the 
4s-component as local potential). Solid lines correspond 
to the all-electron potential, dashed lines to the semilo- 
cal pseudopotentials, and long-dashed lines to their Klein- 
man-Bylander form. Reference energies are marked by solid 
triangles. In the all-electron case the pole in the 4d-channel 
at ~ — 1 hartree is associated with the 3d core state. 
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FIG. 5. Deviations in the excitation energies of the germa- 
nium pseudo atom compared to all-electron results, calculated 
from total energy differences with respect to the groundstate 
configuration. The left panel corresponds to the consistent 
approach using the same XC scheme throughout. The right 
panel refers to the inconsistent approach, using the LDA pseu- 
dopotential but GGA for the XC energy. Solid symbols stand 
for LDA, open symbols for PW GGA values. Squares (□) 
correspond to calculations within linearized CV XC, and tri- 
angles (A) to those within nonlinear CV XC. Results obtained 
within the frozen-core approximation are shown for compari- 
son and marked by circles (o)- Lines are meant to guide to 
the eyes. Note that the underlying excitation energies reach 
up to « 8 eV. 

We have applied criterion (2) employing excited neu- 
tral and (positively) ionized configuations of the spher- 
ical isolated atom. In the case of Ge, e.g., examining a 
4s — ^ 4p electron transfer to mimic orbital hybridization 



upon bond formation, and the first ionization potential. 
In Fig. H we plot the error of the excitation energies for 
consistent calculations within LDA and GGA with re- 
spect to all-electron calculations. 

For linearized CV XC we find the ensuing errors to be 
of the same magnitude albeit of opposite sign as those 
in a frozen-core calculation, where only the all-electron 
valence states are allowed to adjust self-consistently but 
the core charge density is kept fixed as that in the atomic 
groundstate. For nonlinear CV XC the errors of the pseu- 
dopotential calculation approach those of the frozen-core 
calculation. Thus we conclude that the pseudopotential 
related errors are indeed comparably small, absolutely 
and relatively, as those due to neglect of core relax- 
ation, the excitation energy for the transfer (4s^,4p^) — > 
(4s^, 4p^) being about 8 eV. Carrying out the tests in in- 
consistent manner - using GGA XC and LDA pseudopo- 
tentials with linearized CV XC - leads to large devations. 
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FIG. 6. Deviations of the level spacing of the 4s and 4p 
states of the germanium atom with respect to all-electron re- 
sults. See Fig. ^ for a legend. The level spacing is « 8 eV 
for the groundstate configuration, and varies by ~ 0.8 eV. 

Particularly excitation energies are overestimated com- 
pared to full GGA calculations, suggesting that the pseu- 
dopotentials in LDA are more attractive than in GGA, in 
accordance with our findings in Sec. ^ and Employ- 
ing LDA pseudopotentials together with nonlinear CV 
XC we obtain agreement with the GGA frozen-core cal- 
culation however. Considering the eigenvalues we have 
observed an analogous pattern. This is illustrated in 
Fig. H for the deviation of the level spacing eAp — sAs 
with respect to all-electron calculations within LDA and 
GGA. Again a large error occurs in case of the incon- 
sistent GGA calulation using LDA pseudopotentials and 
linearized CV XC, while the consistent approach is accu- 
rate, errors being of the order of a few ten meV compared 
to changes in the level spacing of about 0.8 eV. 

In summary these tests affirm our conclusion reached 
for the bulk systems: Within GGA and LDA the respec- 
tive pseudopotentials possess similar transferability, er- 
rors due to the usual linearization of CV XC being small; 
LDA pseudopotentials are adequate in GGA calculations 
(and vice versa) only if CV XC, acting more repulsively 
in GGA compared to LDA, is incorporated explicitly. 
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TABLE I. Cohesive properties of Na. The first column 
indicates the XC scheme used to generate the pseudopoten- 
tials, the second the one employed for the XC energy of the 
(pseudo) atom and sohd. Bracketed values are based on non- 
linear core- valence XC. We show the lattice constant ao, the 
bulk modulus Bo, and the cohesive energy Et- The latter in- 
cludes a spin correction of the free Na atom of 0.20 eV (LDA) 
and 0.22 eV (GGA). 



potential 




ao (A) 


Bo (GPa) 


Et (cV) 


LDA 


LDA 


3.98 (4.05) 


8.7 (9.1) 


1.28 (1.22) 


LDA 


BP 


3.97 (4.22) 


8.7 (7.3) 


1.06 (0.94) 


BP 


BP 


4.20 (4.22) 


7.3 (7.4) 


0.94 (0.94) 


LDA 


PW 


3.98 (4.21) 


8.7 (7.3) 


1.19 (1.08) 


PW 


PW 


4.24 (4.21) 


7.0 (7.3) 


1.05 (1.08) 


LDA'' 




4.05 


9.2 




PW 




4.22 


7.1 




Experiment'' 


4.23 


6.92 


1.11 



^All-electron data from Reference | 
''Reference 66l 



TABLE II. Cohesive properties of NaCl. Like Table |, and 
using a spin correction of 0.20 eV (LDA) and 0.22 eV (GGA) 
for the free CI atom. 



potential 




ao (A) 


Bo (GPa) 


Et (eV) 


LDA 


LDA 


5.19 (5.43) 


32 (32) 


7.28 (6.99) 


LDA 


BP 


5.27 (5.67) 


27 (22) 


6.59 (6.24) 


BP 


BP 


5.74 (5.68) 


21 (22) 


6.20 (6.19) 


LDA 


PW 


5.28 (5.65) 


28 (23) 


6.85 (6.43) 


PW 


PW 


5.87 (5.66) 


18 (23) 


6.22 (6.40) 


Experiment'' 


5.64 


24.5 


6.51 



'Reference 



TABLE HI. Cohesive properties of hep Mg. Like Table ^ 
The equilibrium c/a ratio was obtained as 1.59 (LDA) and 
1.66 (GGA), the experimental value is 1.62.'' 



potential -Exc 


ao (A) 


Bo (GPa) 


Eb (eV) 


LDA LDA 
BP BP 
PW PW 


3.05 (3.16) 
3.17 (3.18) 
3.20 (3.20) 


39 (37) 
32 (31) 
30 (32) 


2.09 (1.76) 
1.27 (1.22) 
1.42 (1.40) 


Experiment" 


3.21 


35.4 


1.51 



"Reference 

TABLE IV. Cohesive properties of Al. Like Table § includ- 
ing a spin correction of the free Al atom of 0.15 eV (LDA) 
and 0.19 eV (GGA). 



potential 


Exc 


ao (A) 


Bo (GPa) 


Eb (eV) 


LDA 


LDA 


3.97 (3.97) 


83 (85) 


4.09 (4.09) 


LDA 


BP 


3.97 (4.03) 


80 (75) 


3.39 (3.27) 


BP 


BP 


4.05 (4.05) 


75 (75) 


3.26 (3.25) 


LDA 


PW 


3.97 (4.02) 


81 (77) 


3.64 (3.54) 


PW 


PW 


4.05 (4.04) 


79 (79) 


3.52 (3.53) 


LDA" 




3.98 


83.9 




PW" 




4.10 


72.6 




Experiment'' 


4.05 


77.3 


3.39 



"All-electron data from Reference | 
''Reference pfil 



TABLE V. Cohesive properties of diamond. Like Table ^ 
using a spin correction of 1.13 eV (LDA) and 1.26 eV (GGA) 
for the C atom. 



potential 


Exc 


ao (A) 


Bo (GPa) 


Et (eV) 


LDA 


LDA 


3.54 (3.54) 


436 (436) 


8.96 (8.93) 


LDA 


BP 


3.55 (3.58) 


421 (406) 


7.93 (7.57) 


BP 


BP 


3.59 (3.59) 


399 (400) 


7.56 (7.58) 


LDA 


PW 


3.54 (3.58) 


424 (406) 


8.09 (7.92) 


PW 


PW 


3.58 (3.58) 


408 (405) 


7.80 (7.83) 


Experiment" 


3.57 


442 


7.37 


"Reference 








TABLE VI. Cohesive properties of Si. Like Table |, using 


a spin correction of 0.66 eV (LDA) and 0.79 eV (GGA) for 


the Si atom. 










potential 


Exc 


ao (A) 


Bo (GPa) 


Et (eV) 


LDA 


LDA 


5.38 (5.39) 


94 (94) 


5.34 (5.32) 


LDA 


BP 


5.40 (5.46) 


91 (86) 


4.60 (4.47) 


BP 


BP 


5.47 (5.47) 


85 (85) 


4.46 (4.45) 


LDA 


PW 


5.39 (5.45) 


92 (87) 


4.79 (4.66) 


PW 


PW 


5.46 (5.46) 


87 (87) 


4.64 (4.64) 


LDA" 




5.41 


96 


5.28 


BP" 




5.54 


80 




PW" 




5.50 


83 




Experiment'' 


5.43 


98.8 


4.63 



"All-electron data from Reference | 
''Reference 66l 
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TABLE VII. Cohesive properties of Ge. Like Table |, using 
a spin correction of 0.60 eV (LDA) and 0.74 eV (GGA) for 
the Ge atom 



TABLE X. Cohesive properties of bee W. Like Table |, 
using a spin correction of 2.04 eV (LDA) and 2.32 eV (GGA). 



JJULcliLidi J-^XC 


"0 


-Do i^'^r^aj 




LDA LDA 


5 56 (5 60~) 


73 (71) 


4 75 (A 58") 


LDA BP 


5.59 (5.74) 


67 (57) 


3.96 (3.66) 


BP BP 


5.73 (5.76) 


59 (56) 


3.70 (3.66) 


LDA PW 


5.58 (5.74) 


69 (59) 


4.14 (3.82) 


PW PW 


5.74 (5.74) 


58 (58) 


3.82 (3.82) 


LDA'' 


5.63 


78 


4.54 


gpa 


5.76 


60 




PW 


5.75 


61 




Experiment'' 


5.66 


76.8 


3.85 


^All-electron data from Reference 


28[ 




Reference |66[ 








TABLE VIII. Cohesive properties of GaAs. 


Like Table El 


using a spin correction of 0.15 eV (LDA) and 0.18 eV (GGA) 


for the Ga, and 1.41 eV (LDA) and 1.67 eV (GGA) for the 


As atom. 








potential -Exc 


ao (A) 


Bo (GPa) 


Et (eV) 


LDA LDA 


5.50 (5.57) 


79 (75) 


8.68 (8.15) 


LDA BP 


5.52 (5.72) 


74 (62) 


7.19 (6.33) 


BP BP 


5.68 (5.72) 


64 (62) 


6.52 (6.33) 


LDA PW 


5.51 (5.70) 


76 (64) 


7.51 (6.63) 


PW PW 


5.69 (5.71) 


63 (64) 


6.74 (6.63) 


LDA'' 


5.62 


74 


7.99 


gpa 


5.76 


60 




PW 


5.74 


65 




Experiment'' 


5.65 


74.8 


6.52 


''All-electron data 


from Reference 


2| 




''Reference [jil 








TABLE IX. Cohesive properties of fee Cu. 


Like Table ffl 


using a spin correction of 0.20 eV (LDA) and 0.25 eV (GGA). 


potential -Exc 


ao (A) 


Bo (GPa) 


Et (eV) 


LDA LDA 


3.55 (3.56) 


172 (172) 


4.31 (4.24) 


LDA BP 


3.68 (3.70) 


124 (122) 


3.09 (3.06) 


BP BP 


3.67 (3.68) 


130 (131) 


3.22 (3.23) 


LDA PW 


3.67 (3.69) 


127 (123) 


3.23 (3.20) 


PW PW 


3.67 (3.67) 


134 (132) 


3.38 (3.38) 


LDA" 


3.52 


192 


4.29 


BP" 


3.62 


151 


3.12 


PW" 


3.62 


151 


3.30 


Experiment'' 


3.60 


138 


3.50 



potential 


-Bxc 


ao (A) 


Bo (GPa) 


Et (eV) 


LDA 


T T— \ A 

LDA 


3.14 (3.16) 


324 (331) 


10.75 (10.24) 


LDA 


BP 


3.15 (3.21) 


308 (306) 


9.26 (8.52) 


BP 


BP 


3.20 (3.21) 


299 (308) 


8.73 (8.54) 


LDA 


PW 


3.14 (3.20) 


313 (308) 


9.63 (8.87) 


PW 


PW 


3.22 (3.21) 


298 (310) 


8.88 (8.87) 


LDA" 




3.14 


337 




PW" 




3.19 


307 




Experiment'' 


3.16 


310 


8.90 



"All-electron values from Ref. 
''Reference pfil 



TABLE XI. Change of the cohesive energy, SEt, due to the replacem 
energy change per atom as obtained from self-consistent calculations, SI: 
The latter is decomposed into its consituent terms arising from the d 
pseudopotentials, SEy. The values without and with brackets corresponc 
nonlinear core- valence XC respectively. No spin-corrections were applied ti 
at their experimental values in all calculations, the ensuing error of the t 
equilibrium structure staying below ~50 meV per atom. Symbols are exf 



r TT^solid 
O-Etot 



(eV) 



Na 


SEtot 


0.32 (-3.21) 




SEx_c + SEv 
SExc 


0.32 (-3.20) 
-0.01 (-3.20) 
0.33 (-0.00) 


Diamond 


SEtot 


0.31 (-1.42) 




6Ex_c + SEy 
SExc 
SEv 


0.32 (-1.40) 
-0.46 (-1.37) 
0.78 (-0.03) 


Si 


SEtot 


0.60 (-0.97) 




SEx_c + SEy 
SExc 
SEy 


0.61 (-0.96) 
-0.27 (-1.09) 
0.88 ( 0.13) 


Ge 


SEtot 


1.50 (-5.50) 




SExc + SEy 
SExc 
SEy 


1.50 (-5.48) 
-0.24 (-5.49) 
1.74 ( 0.01) 


Cu 


SEtot 


1.63 (-1.99) 




SExc + SEy 
SExc 
SEy 


1.64 (-1.99) 
-5.54 (-7.48) 
7.19 ( 5.50) 



" All-electron values from Ref. 
''Reference 66l 



(ao, Bo), and Ref. y (Eb). 
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